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I. INTRODUCTION 

Weakly interacting ultracold Bose gases are of great 
interest because they allow the 'pure' study of quantum 
phenomena on a macroscopic scale. Such systems are 
typically engineered in harmonic traps with different ge- 
ometric configurations, whose dimensionality plays a cru- 
cial role in determining the properties of these gases. In 
three-dimensional harmonic traps, the system undergoes 
a phase transition which leads to the appearance of coher- 
ence across the entire sample, as demonstrated in (THl]- 
By setting the trap frequency in one direction to be much 
larger than the others, the effective system dynamics is 
reduced to two dimensions 5 8] , with many interesting 
phenomena occurring, such as the Berezinskii-Kosterlitz- 
Thouless transition [9Nl6j| . Increasing the trap frequency 
in a further direction allows the realization of highly elon- 
gated traps, where the interesting physics occurs in the 
axial direction [|, [l7l423j . 

In such a regime two characteristic temperatures be- 
come relevant, associated with the onset of phase (T^) 
and density (Td) fluctuations [24|. For temperatures 
< T < Td, density fluctuations tend to be sup- 
pressed, and the system reduces to a condensate with 
fluctuating phase (quasicondensate) [25j]. In such a ge- 
ometry low-energy thermal excitations of the axial modes 
play a crucial role, as they tend to destroy the coher- 
ence in the sample |24l - [40j ; such excitations may have 
wavelengths greater than the transverse extent of the sys- 
tem j_therefore acq uiring a one-dimensional (ID) charac- 
ter |19l-[23. [3a - t49| . An accurate analysis of the coher- 
ence properties in such systems is therefore necessary for 
potential a ppli cations, such as matter-wave interferome- 
try (HI . [50l455| | , atom chips [56| and atom lasers 0, [57r - 
|60| . Although experiments can nowadays be engineered 
to produce gases which are both weakly interacting and 
practically ID (/j, fc^T < hjJ±, where (i is the chemical 
potential, ksT the thermal ener gy a nd Hlo± the trans- 
verse excitation energy) (23l . l4lT 49] . the early experi- 
ments performed did not satisfy these conditions so well, 



and the system was instead in the 1D-3D crossover regime 

In this work we propose and implement a modified 
stochastic model which enables us to perform a suc- 
cessful ab initio description in the 'intermediate' regime 
/i, ksT ~ few huj± where quasi-condensate physics dom- 
inates, but transverse effects still need to be appropri- 
ately accounted for. More specifically, we analyze two 
early experiments [13, EO] investigating the phase co- 
herence properties of weakly interacting, quasi- ID Bose 
gases. These experiments measured the temperature de- 
pendence of the coherence length L c , over which the gas 
maintains an appreciable coherence, both in the 'strong' 
(6 < T/T < 28) Q and ' wcak ' (0-8 < T/T$ < 8) 
[4(| phase fluctuation regimes. Although it is theoreti- 
cally anticipated (for a homogeneous gas [6l[), that the 
coherence length, scaled to the experimental half-length 
of the system L, should yield a universal curve when 
plotted against the reduced temperature T /T^, it would 
not actually be appropriate to incorporate their reported 
data into a single graph [62j spanning the entire regime 
0.8 < T/T^ < 28, because the two experiments used dif- 
ferent techniques to extract the coherence. 

The main achievements of this work are as follows: 
(i) we show that the measurements reported in both 
experiments are indeed consistent with such a unified 
curve; this is achieved by a completely ab initio analysis 
of the experimental data, using only the quoted exper- 
imental values for atom number, trap frequencies and 
temperature; (ii) our analysis, which is based on the 
quasi-ld stochastic Gross-Pitaevskii equation [g3|, ap- 
propriately modified here to extend its validity to the 
regime /i ~ few huj± reveals excellent agreement in the 
'strong' phase fluctuation regime; (hi) our attempt to 
more closely mimick the experimental procedure of (40j 
- whose findings have not been adequately interpreted 
to date - also yields good overall agreement (within er- 
ror bars); in order to further improve this we implement 
an alternative approach of extracting a phase coherence 
temperature based on direct analysis of the phase distri- 
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butions from different realisations of our stochastic treat- 
ment. 



II. METHODOLOGY 

Our treatment is based on a modified one-dimensional 
form [63| of the stochastic Gross-Pitaevskii equation 
(SGPE) The modification proposed in - 

which was found to be essential to accurately simultane- 
ously reproduce both in situ density profiles and density 
fluctuations of recent quasi-one-dimcnsional experiments 
El IM El - takes the form: 
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+Huj x [y/1 + 4\ip\ 2 a s - 1] - V»C*, *) + ^0, *)■ (1) 

Here tp describes the relevant (highly populated) low- 
energy modes of the system, V(z) = mLo 2 z 2 /2 is the 
axial trapping potential and rj is a complex Gaus- 
sian noise term, with correlations (rj* {z,t)rj{z' \t')) — 
2h'f(z,t)kBTS(z — z')S(t — t'), where 7(2, t) denotes the 
dissipation. (For further details on numerical implemen- 
tations and related treatments see Refs. f28l Ir33l - |8G ] ) . The 
condition fi < few fku± can lead to a swelling of the con- 
densate in the transverse direction, relative to the true ID 
transverse ground state; this quasi-lD effect is due to re- 
pulsive interactions and is reproduced by using the mod- 
ified nonlinear term huj_ [yl + 4|^'| 2 as — l] (where a s is 
the s-wave scattering length), which reduces to the ID 



result in the limit 4|^>| a s <C f . 



For the ordi nary Gross- 



Pitaevskii equation, this was shown in Refs. 
and its validity was verified experimentally in [48f . 

In order to match to experimental atom numbers, it is 
crucial to also include in our treatment the contribution 
(n±) to the linear density profile of transverse thermal 
atoms with energy greater than fku± . This is done via 



n(z;n,T) = (\ip(z; fi, T)\ 2 ) + n±(z;fx,T) 



(2) 



where (• • • ) denotes ensemble averaging, obtained by av- 
eraging over many different realizations of the noise, and 
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n^z-^T) = —^( J + l)g 1/2 



(3) 



where ffi/2(---) is the polylogarithm (or Bosc function) of 
order 1/2, and XdB = h/ \J2nmkBT is the thermal de 
Broglie wavelength J85|. While in the limit \x < tvu)\_ 
studied in our previous work (63j . it is sufficiently accu- 
rate to use fij{z) = fi — V(z) — jhu± in Eq. §5§ (see also 
[UllisI); the effect of mean- field potential experienced by 
the transverse thermal atoms should also be taken into 
account in the regime [i > characteristic of the ex- 
periments of Refs. [22], |4(| studied here. For this reason 
we use here the modified expression: 

Hj(z) = ft- V(z) - jhu ± - 2. 9 ((|Vf) + n±) (4) 



where g = 2a s Hui± is the effective interaction strength. 

In the following sections we systematically use 
Eqs. (HJ-BJ to match to the experimental atom number 
(noting also that the profiles generated by Eq. ([2]) corre- 
spond to transversely integrated density profiles, typical 
of ultracold gas experiments), before performing an anal- 
ysis of the coherence properties. 

For both experiments analysed here, the only experi- 
mental inputs to our theory are the atomic species used 
( 87 Rb in both cases), the trapping configuration, the 
temperature and atom number corresponding to each re- 
ported experimental data point (87j . In our ab initio anal- 
ysis we fix experimental trap configuration and temper- 
ature, and vary /1 until the desired number of atoms is 
reached in our simulations. Depending on the available 
experimental data, we either match to the total or the 
quasicondensate experimental atom number, as explained 
in each case. Representative error bars are also calculated 
for our simulations, based on variations in both number 
and temperature, in close analogy to the experimental 
analysis. 

In order to discuss both experiments in a coordinated 
way, and demonstrate their consistency, we consider the 
dependence of the coherence length on temperature, with 
both parameters appropriately scaled to show the emer- 
gence of universal physics. In performing this analysis, 
we also address experiment-specific features, meaning the 
precise definitions of the coherence length and phase fluc- 
tuation temperature differ slightly within our subsequent 
analysis, however this difference is clearly labelled within 
figures. 



III. COMPARISON TO THE EXPERIMENT OF 
RICHARD ET AL. 

The first experiment we consider is that of Richard 
et al. [23], in which a very elongated harmonic trap 
(lu±/2tt = 760 Hz and uj z /2tt = 5 Hz) was used to gen- 
erate quasicondensates of atom numbers in the range 
0.25 x 10 5 — 0.65 x 10 5 and at temperatures 90 nK 
< T < 350 nK. These parameters enabled the 'strong' 
phase fluctuation regime 6 < T/T^ < 28 to be probed. 

In order to investigate the phase coherence properties 
of the gas, the axial momentum distribution was mea- 
sured by means of Bragg spectroscopy and found to have 
a Lorentzian shape, consistent with an exponential de- 
cay of the correlation function in space anticipated in the 
regime of large phase fluctuations [2l|, [22|, Hi| HH . This 
approach enabled the coherence length L c to be extracted 
from the half width at half maximum of the momentum 
profile, and its temperature dependence to be measured. 
Using the original data relevant for Ref. [22| (see also [39| 
and |88(). we can plot this in universal form by scaling 
the coherence length L c to the half-length of the con- 
densate L, and plotting against the reduced temperature 
T/T,),. Here we use the general definition of in terms 
of the ID axial quasicondensate peak density, n qc (0), i.e. 
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Figure 1. (Color online) Coherence length L c scaled to the 
half length of the quasicondensate L vs. reduced tempera- 
ture T/T^. Comparison of experimental data [22l [39(] span- 
ning the range T /T^ > 6 (hollow red triangles) to the SGPE 
model (upward filled black triangles). Additional numerical 
results generated with the SGPE model for lower tempera- 
tures T/T$ < 6 are also shown (downward filled black tri- 
angles), together with the best fit of the SGPE results (dot- 
dashed black line). The horizontal experimental error bar 
is taken as 10% [2Iq|. based on typical experimental uncer- 
tainties, while the corresponding error in the simulated point 
arises from a typical 15% variation in the quasicondensate 
atom number and 5% variation in temperature [22|; in both 
cases the vertical error bars fall within the symbol size. Tem- 
peratures have been scaled to T0[n qc (O)] = fi 2 %(0)/mfcsl, 
where n qc (0) is the peak quasicondensate density; in the 
SGPE analysis, the spatial extent of the quasicondensate is 
obtained ab initio as the temperature dependent Thomas- 
Fermi radius Rtf(T) 



T0[n qc (O)] = h 2 n qc (0)/mk B L , 



(5) 



which is relevant for Refs. |22|, [39| ■ (See discussion fol- 
lowing Eq. (7) for a closely related expression valid in 
the Thomas-Fermi regime and used in the analysis of 
Ref. [H.) 

In order to simulate the experiment, we compute 
the first order correlation function in position space 
gW(-z/2,z/2) = (V>*(-z/2)^(z/2)} (normalized to the 
averaged central density (|^(0)| 2 )), finding an exponen- 
tial behaviour against distance z as in the experiment 
(since T ^> here) . We then compute the Fourier trans- 
form = Flg^] and extract the coherence length L c 
by measuring Afc, the half-width at half maximum of 
CW. We also use our simulated results to find T^,[n qc (0)], 
with the required inputs being the quasicondensate peak 
density and spatial extent of the gas, each of which we 
obtain ab initio. Specifically, the quasicondensate density 
is extracted as [H Ill-Ill : 



n qc (z) = vlWl^f 



(6) 



and its peak value is used in the expression for T^[n qc (0)]. 
Moreover, the half-length of the quasicondensate in 
our simulations is given by a temperature dependent 
Thomas-Fermi radius Rtf{T), extracted from Eq. ([1]) 
analogously to the procedure used in the modified Popov 
theory (see HI HI ill and (H). 

Our ab initio SGPE results presented in Figure [T] (up- 
ward filled black triangles) show excellent agreement to 
the experimental results (hollow red triangles) in this 
'strong' phase fluctuation regime within their respective 
characteristic error bars (see caption for details). Note 
that the coherence length is smaller than the quasicon- 
densate extent (L c /L < 1) in the temperature range 
probed here, illustrating the fundamental role of phase 
fluctuations in such an elongated geometry. As T /T^ ap- 
proaches zero, we expect the coherence length to increase: 
to show this, we have generated a further set of numeri- 
cal points for T/T^ < 6 (downward filled black triangles), 
which indeed confirm this picture, and the 'universal na- 
ture' of such a scaled diagram. 

In the following section, we turn to the investigation of 
the experiment of Hugbart et al. [40j | , in which the regime 
T/Tfj, < 8 was probed using different methods. 



IV. COMPARISON TO THE EXPERIMENT OF 
HUGBART ET AL. gg] 

The Bragg spectroscopy method used in the previous 
experiment limited accurate investigations of the coher- 
ence properties of the gas to T/T^ > 6, as at lower tem- 
peratures the width of the momentum distribution of the 
gas was no longer easily resolved |39j . In the experiment 
described in [4 0| . an alternative interferometry technique 
(see also Ref. [2lj] ) was used to measure the spatial cor- 
relation function in the regime 0.8 < T/T$ < 8. 

In particular, after the condensate was released, two 
Bragg pulses were applied, playing the role of matter- 
wave beam splitters. The contrast of the resulting in- 
terference fringes was then obtained from the modulus of 
the Fourier transform of the interference pattern, and the 
coherence length was extracted from the decrease of the 
contrast as a function of the distance between the two 
interfering condensates. The experimental data (hollow 
red circles of Fig. [5]), show that, in the regime T ~ T$, 
the coherence extends over more than half of the system 
size. 

In order to access such low values ofT/T^, it was tech- 
nically easier to use slightly less elongated traps than in 
[22^ . In fact, two different trap configurations were used 
(first: uj±/2ir = 395 Hz and u> z /2ir = 8.67 Hz; second: 
co±/2ir = 655 Hz and oj z /2it = 6.55 Hz); the data ob- 
tained with the second trap were subdivided into two 
different blocks, characterised by two different values of 
the evaporation parameter, which proved necessary in or- 
der to vary T/T^, while keeping the condensed fraction 
fairly constant within each data block. The total atom 
numbers measured in this experiment were found to lie 
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Figure 2. (Color online) Scaled coherence length L c /L vs. re- 
duced temperature T/T<j,. Comparison of experimental data 
[4Cj | (hollow red circles) to predictions of the SGPE model 
(filled black squares). The reported experimental best fit 
(dashed red line) and the theory used in [4(J to interpret the 
experimental data (solid line) are also shown, alongside the 
SGPE fit from Fig. Q] (dot-dashed black line). Typical error 
bars for corresponding data points in both experiment and 
SGPE model are indicated, with the SGPE ones arising from 
a characteristic 20% variation in total atom number and an 
additional 5% variation in temperature; (note that the SGPE 
error bars for the point with T/T$ ~ 1.2 lie within the point 
size and are not visible). The SGPE results are scaled to 
2^>[n q c(0)] for consistency with the discussion of Fig. [1] 



in the range 0.8 x 10 5 — 3 x 10 5 (corresponding quasicon- 
densate [93| atom numbers: 0.5 x 10 5 — 2.5 x 10 5 ), within 
a temperature region of 100 — 230 nK |87j . 

In our simulations, rather than reproducing the ex- 
perimental procedure, for which non-equilibrium expan- 
sion dynamics would need to be accounted for, we ini- 
tially instead extract the corresponding in situ coherence 
length by adopting the same methodology described in 
Section UTTl This is done here in order to firstly explore 
whether our results for this new system also lie on the 
same 'universal' curve as those of Ref. [22[ (dot-dashed 
black curve of Fig. Q]) when analysed in an identical man- 
ner. 

We indeed find that the new numerical results gen- 
erated in this regime T ~ fewT^ (filled black squares of 
Fig- HI) lie on the same (dot-dashed black) curve provided 
by the fit in Fig. [TJ this demonstrates the universal char- 
acter of the coherence properties, as numerical results 
obtained with different sets of trap configurations, tem- 
peratures and atom numbers yield the same behaviour. 
As the analysis of Ref. [40| reported error bars for two in- 
dicative data points, we repeat a similar analysis within 
the SGPE theory - in which case it is based on a 20% 
variation in total atom number [94| . and a typical 5% 
variation in temperature (95l . l96j . 



However, while our methodology led to a unified theo- 
retical graph over both the 'weak' and 'strong' phase fluc- 
tuation regime, which also provides excellent agreement 
with the experiment [22j for T > BT^, the correspond- 
ing SGPE results for the experiment of Ref. [4(j differ 
from the experimentally-reported ones. The reason for 
this discrepancy is twofold: firstly, in [4(| it is an 'effec- 
tive' correlation function that is evaluated, which leads 
to a different definition of the coherence length [62[ (see 
Section HV Aft ; secondly, the theoretical values of T/T^ 
obtained in our treatment do not span precisely the same 
range as in the experiment (see Section IIVBI) . 

With regard to the latter point, we note that the anal- 
ysis of Ref. [Io| was actually based on a slightly modified 
definition for compared to that of Eq. (5); in par- 
ticular, in [40j . was defined in terms of the number 
of quasicondensate [§3[ atoms A qc (and not on the peak 
quasicondensate density), via the expression (97j 



7>[A qc ] = 15H 2 N qc /l6mk B L 2 



(7) 



This is in fact a simplified form of Eq. (5), valid for 3D 
condensates [98| (i.e. condensates where the transverse 
density profiles can be well approximated by a Thomas- 
Fcrmi profile). 

In the following sections, we further investigate poten- 
tial sources of discrepancy between our simulations and 
the experimental results of Ref. |4(j, focusing our anal- 
ysis on a modified effective correlation function defined 
below, as relevant for this experiment. Throughout our 
analysis, the coherence length is always referred to as 
L c , with the method of its extraction clearly identified in 
each figure. 



A. Effective correlation function 

The main source of the observed discrepancy should 
be related to the experimental measurement of an 'effec- 
tive' correlation function (see Eq. (9) in |4(J), instead of 
the correlation C' 1 ' measured in Richard et al. [22[ and 
discussed thus far. In the experiment of Ref. [40| . this ef- 
fective correlation was introduced in order to cancel the 
random phase caused by the shot to shot fluctuations of 
the global position of the contrast fringes; it was found 
that taking the absolute value of the Fourier transform of 
the fringe pattern before averaging achieved this aim (for 
a more detailed explanation see Section 4.2 in Ref. j40|). 
but modified the coherence length relative to that of C^- 1 ' . 
For this reason, we should not expect the measurements 
from the two experiments to lie on the same curve, as 
they measure two different quantities [62| . 

In analogy to the experimental method used to ex- 
tract the coherence length, we implement this feature 
by taking the absolute value of the Fourier transform of 
g^ 1 ' from each individual run, before averaging over the 
different realizations of the noise; our 'effective' corre- 
lation function, which we here call C^' mod \ takes the 
form C( 1 > mod ) = (\F[g^(-z/2,z/2)}\). The 'effective' 
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correlation function is found to have similar behaviour 
to but it decays faster (in the region of interest) in 
momentum space, thus resulting in larger values of the 
coherence length. This is shown in Fig. [3J which com- 
pares the theoretical results to the experimental mea- 
surements using C^ 1 - 1 (brown filled squares) and C'( 1 ' mod - ) 
(blue filled circles). Although we consciously do not ex- 
actly reproduce the experimental measurement sequence, 
the results obtained from the SGPE analysis of the 'effec- 
tive' correlation function clearly show a very similar trend 
to the experimental findings over the probed regime; in 
particular, they tend to lie on the reported line of best 
fit of the experimental data (dashed red line). We note 
that our calculation of C'( 1 ' lnod ) leads to a much improved 
agreement with the experimental data than the original 
theoretical analysis reported in poj (solid red line) which 
largely overestimates the amount of coherence in the sys- 
tem. 




Figure 3. (Color online) Scaled coherence length L c /L vs. 
reduced temperature T/Tj,. Comparison of numerical SGPE 
findings extracted via C (1 ' mod) (filled blue circles) and via 
C' 1 ' (filled brown squares) against experimental data points 
of Ref. [4(| (hollow red circles). The theory used in [40T | is 
also shown (solid line) together with the experimental best 
fit (dashed red line). SGPE numerical points are scaled to 
T^>[iVqc], and respective error bars account for 20% variation in 
total atom number and a further 5% variation in temperature. 
Inset: comparison of SGPE results, in which temperatures are 
scaled to the experimentally quoted values of . SGPE results 
shown in the inset are obtained by matching either the total 
number of atoms (blue filled circles, main figure and inset), or 
the quasicondensate atom number (green crosses, inset only). 



Although the experimental results and our simulated 
points (extracted via C (1 ' mo ^) demonstrate very good 
agreement when accounting for their respective error 
bars, the experimental data still appear to systemati- 
cally extend to slightly larger values of T/T^, as visible 



in Fig. [3j This could be attributed either to a systematic 
shift in the experimental determination of T (e.g. aris- 
ing in expansion imaging), which however increases with 
increasing T '/X^ ratio, or to the method by which the 
inputs to Tfi are extracted in the analysis. In the rest 
of the paper we assume that this shift arises solely from 
the latter and attempt to further improve on the spanned 
range of T/T . 

B. Identification of 

A given set of trap frequencies, temperatures and to- 
tal atom numbers should fix the characteristic temper- 
ature T,), to a particular value which is the same be- 
tween theory and experiment; for this reason, in order 
to attempt to span precisely the same range of T/T$ as 
in the experiments, we plot in the inset of Fig. [3J our 
theoretically-obtained values for the reduced coherence 
length L c /L versus T/T$ using here the experimentally- 
obtained values for T^. This still reveals good agree- 
ment between our theoretical findings (blue filled circles) 
and the experimentally-extracted ones (hollow red cir- 
cles), within the reported experimental error bars. 

We have also performed a separate analysis, based 
again on the SGPE but in which we instead match our nu- 
merically extracted quasicondensate number to the cor- 
responding experimentally-extracted 'condensate' num- 
ber |87j . and measure the coherence length via 
These closely related measurements are shown in the in- 
set of Fig. [3J with green crosses; we find good overall 
agreement both with the experimental findings (hollow 
red circles) and with the previous analysis, based on 
matching the experimental total atom number instead 
(filled blue circles). 

In order to further resolve the remaining discrepancy 
between theory and experiment in the regimes of T ' /T$ 
analysed, in the next section we discuss an alternative 
method to extract from the SGPE simulations. 

1. T,p extracted from the phase distribution 

In this section we investigate an alternative method of 
reproducing the experimental results of Ref. (ifjj from an 
SGPE analysis. Our approach is motivated from foot- 
note 47 of [4(|, stating that can be obtained from 
the relation L^/L = T^/T [l7| with identified as the 
characteristic separation over which the phase fluctuates 
by 1 radian at the trap centre. 

In order to extract values for according to the rela- 
tion described above, we take here the approach of sys- 
tematically analysing the phase distributions of the en- 
semble of stochastic fields ij; at several distances from the 
trap centre. Within the Thomas-Fermi radius, we find 
these distributions to be well fitted with Gaussian func- 
tions, whose standard deviation increases with increasing 
distance from the trap centre due to the enhanced role of 
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thermal fluctuations (see Fig. [5] in Appendix [A] and re- 
lated work in |47|)- In this analysis, can be identified 
as the distance from the centre where the standard devi- 
ation of the Gaussian fit to the phase distribution reaches 
a particular value. Our analysis indicates that when this 
(free) parameter takes the value of a = 0.65, optimum 
agreement with the experimental findings is obtained re- 
garding the spanned range of T/T^. 
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Figure 4. (Color online) Scaled coherence length vs. re- 
duced temperature. Experimental data [4(| (hollow red cir- 
cles) against numerical SGPE results extracted via 

c (l, m od) 

(blue stars); in the SGPE analysis temperatures are scaled 
to extracted from the phase distribution (see Appendix 
|A")) . The theory used in [ZFJl is also shown (solid line) together 
with the reported experimental best fit (dashed red line). Er- 
ror bars account for 20% variation in total atom number (solid 
blue line) and 8% variation in the standard deviation (dashed 
blue line) (see Appendix [XJ. 

The simulated results, with temperatures scaled to 
extracted as described above (blue stars), are shown in 
Fig. |U alongside the experimental findings (hollow red 
circles). We also include here error bars for the previ- 
ously considered indicative data points (as in Ref. (40j|). 
In obtaining these we consider here two sources of er- 
ror: firstly, as in previous figures, a variation in total 
atom number of 20% 94|; secondly, as the horizontal 



axis is also sensitive to the precise value of a chosen, the 
sensitivity to a variation in a of ±0.05 is shown by the 
additional dashed portions of the error bars. 



model, the stochastic Gross-Pitaevskii equation (SGPE), 
with novel and self-consistent modifications in order to 
take account of quasi-one-dimensional effects. Compar- 
ing to the experiment of Richard et al., Ref. [22 |. which 
considered the regime of relatively strong phase fluctua- 
tions T/Tfi > 6, we found excellent agreement between 
experiment and the SGPE theory, based solely on using 
the experimental parameters (i.e. trap configurations, 
atom number, temperatures) as inputs for the theory. 

We further compared the SGPE theory to the findings 
of Hugbart et al., Ref. [4(| for the opposite low temper- 
ature regime T < few T^, which is more challenging to 
probe in experiments. While the SGPE analysis gave a 
temperature dependence of the coherence length in quan- 
titative agreement with the experimental trend (within 
error bars), undertaking a point-by-point analysis of the 
experimental data was found to span a slightly narrower 
range of T/T^ than the experimentally-reported curve, 
indicating a systematic deviation. We argued that this 
discrepancy may arise as a result of the different identi- 
fications of between theory and experiment, possibly 
due to the different means of processing the 'raw' experi- 
mental data and stochastic numerical results, e.g. due to 
differences in extracting the quasicondensate atom num- 
ber which then feeds into the expression for T^. 

We have partially examined this issue by using in- 
stead a phase sensitive means of extracting from 
characterisation of the ensemble phase distribution in 
SGPE simulations. In particular, motivated by foot- 
note 47 of Ref. [io| , we identified T$ through the relation 
T$ = T(L^j 'Rtf(T)) where was chosen as the char- 
acteristic separation from the trap centre at which the 
phase distribution is fitted by a Gaussian with a partic- 
ular value for the standard deviation. In our treatment 
however this value is a free parameter chosen here so as 
to match the experimental range of T/T^. 

Finally, as the measurements performed in [4(| rely on 
an intcrferometric technique, this opens up a way to indi- 
rectly 'extract' the Penrosc-Onsager mode in a quasi-one- 
dimcnsional experiment solely from knowledge of density 
and the two lowest equal-time correlation functions, with 
the Penrose-Onsager condensate mode density emerging 
as the phase- and density-fluctuation suppressed part 
of the density via iipo(z) ~ n(z) \Jl — g^{z) </W(0, z) 
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Appendix A: Analysis of the phase distribution 

The numerical procedure implementing the SGPE 
model makes use of the Condor system [991 ] , which is a 
mechanism to run parallel simulations, on different ma- 
chines. Using this method we generate a large number of 
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Figure 5. Distribution of the phase of a*ijj(z)/\a c \ (where a c 
is the amplitude of the Penrose-Onsager mode), for a chosen 
temperature T ~ (left point with error bars in Figs. I2I4II . 
at different locations from the trap centre. Shown are the 
distributions at the condensate edge (left plot), at the point 
2 = L,), (central plot) where the distribution is fitted by a 
Gaussian (illustrated by the red solid curve) of standard de- 
viation a — 0.65 (which for the particular numerical point 
considered occurs at z ~ 0.5 7?tf(T)) and at the trap cen- 
tre z — (right plot). The distributions are centered and 
normalised, and the phase <j> is scaled to n. 



stochastic realizations of the wavefunction tp (typically 
~ 1000 in a very short time); in this way we can then 
extract physical observables by performing the average 
over the ensemble of stochastic trajectories. In each nu- 
merical realization, which is somewhat analogous to an 
individual experimental run, the phase of the wavefunc- 
tion tp takes on a certain value at each position, leading 
to a distribution across the ensemble; it is then inter- 
esting to investigate how the distribution of the several 
realizations of the phase changes with the distance from 
the trap centre, at a definite value of T/T$. 

In Fig. [5] we report phase histograms of the stochastic 
field ip(z), locked to the p hase of the Penrose-Onsager 
(condensate) mode <f> c [ZJ, Il00| . at three values of the 
distance from the trap centre, for the experimental data 
point with T « (left point with error bar in Figs. [5J- 
[4}. To be more specific, we plot the phase of a*ip(z)/\a c \, 
where a c is the amplitude of the Penrose-Onsager mode, 
given by a c = A z ^2 7 . <j)^(zi)il)(Zi), where Az is the grid 
spacing (see [73, Il00| for further details and implemen- 
tation). 

The broadness of the generated distributions is an in- 
dication of the amount of coherence at a specific spa- 
tial point in the system: we expect the distribution to 
become broader with increasing distance from the trap 
centre, and be almost flat at the edge, where the system 
becomes purely thermal (Fig. [5]). We define L$ to be 
the distance at which the phase histogram is well fitted 
by a Gaussian with a particular value of the standard 
deviation a. Following such a procedure, we can then 
define a phase fluctuation tem per ature from the rela- 
tion T w T(L /i?TF(T)) [Hill]. The standard devia- 
tion is here a free parameter, and we find that the value 
of a = 0.65 ± 0.05 provides very good matching to the 
experimental data, as shown in Fig. U (with correspond- 
ing range z = ~ (0.15 — 0.6) Rtf{T) for this chosen 
value of a). 
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